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1. Introduction 

Many interesting statistical-mechanical systems can be described by 1+1-dimensional 
models with boundaries. These include the Kondo problem and the Anderson model, tun- 
neling in quantum wires, dissipative quantum mechanics, and the Callan-Rubakov effect. 
At and near a critical point, they can be studied using boundary conformal field theory 
J]]]. Moreover, many of these models are integrable and many properties can be derived 
exactly even away from critical points |§. 

There are several approaches to integrable models, with and without a boundary. One 
is explicit diagonalization of the (continuum or suitably discretized) Hamiltonian using the 
Bethe ansatz. Another bypasses this diagonalization by using maximally the constraints 
resulting from integr ability: after making some guesses for the particle content one can 
find physical quantities like the exact S matrix for particles scattering among themselves 
and off the boundary. A crucial ingredient in both cases is the Yang-Baxter equation, 
including the boundary part ||. In the second approach this equation only allows one to 
find ratios of S matrix elements and not their overall prefactor. In the bulk such prefactors 
are controlled by unitarity and crossing symmetry. For the boundary, constraints of the 
same type were found in , enabling the determination of S matrix elements up to the 
usual CDD-type ambiguities. This was done for the Ising model in a boundary magnetic 
field ||, the sine-Gordon model with a boundary potential affine Toda field theories 
|0,H and the 0(n) sigma models with free or fixed boundary conditions ||. 

The purpose of this paper is to use the first approach to confirm the second. We extract 
boundary S matrices from lattice regularizations of integrable quantum field theories, 



generalizing a well-known method in the bulk [10,11.12]. Our explicit example is the 
inhomogeneous XXZ model in a boundary magnetic field, which is expected to provide an 
integrable regularization of the sine-Gordon model with fixed boundary conditions. Our 
results confirm results of [|J, in particular their cross- unitarity relation, and remove the 
CDD ambiguity. The interest of the computation lies mainly in the fact that boundary 
effects are subleading, and a number of subtleties arise in their analysis. 

We derive the boundary S matrix elements in the usual manner. We start with 
the explicit diagonalization of the appropriate lattice model with integrable boundary 
conditions. This results in the "bare" Bethe ansatz equations, which relate the level 
densities of the bare particles to the actual filled densities of these particles. Using these 
equations, we find the ground state of the system, fill the Fermi sea, and identify physical 
excitations. This results in the "physical" Bethe ansatz equations, which involve the 
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densities of the actual quasiparticles of the system. On the other hand, physical Bethe- 
ansatz equations can also be derived directly from the conjectured S matrix, without 
reference to a lattice. We perform such a computation for the S matrix of [f|||. We 
compare the two results, and find that they are indeed the same. 

The lattice model of interest is the inhomogeneous 6-vertex model [|l3l , |l4[l on an open 



strip, with integrable boundary conditions as in |15[]. We take identical boundary conditions 
on both sides of the system. In the hamiltonian and homogeneous limit, this model reduces 
to the XXZ model with boundary magnetic fields h JTB1 



H 



1 
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In the inhomogeneous 6-vertex model, one gives an alternating imaginary part ±iA to the 



spectral parameter on alternating vertices [[13| , p~4|| . The scaling limit is given by taking 
A — > oo, N — > oo, and the lattice spacing A — > 0, such that L = NA remains finite. In 
the bulk, this provides a regularization of the sine-Gordon model with Lagrangian 



L SG = / dx [\{d4>) 2 + n 2 cos (3 sa(f\ 
Jo 



(2) 



where the mass fi oc exp(— const A), and /3| G = 8(^ — 7) for the antiferromagnetic case 
(e = —1), while = 87 for the ferromagnetic one (e = 1). The scaling limit of the 
gapless hamiltonian (p is given by the ultraviolet limit of (0). 

Free boundary conditions in the lattice model correspond to fixed ones in the sine- 
Gordon model. This can easily be seen by recalling that <p is dual to the arrows of the vertex 
model; since in the approach of the time direction is across the diagonal of a vertex, 

the field at the boundary is necessarily constant. Thus we have <p(0) = 4>{L) = </>o, where 
(f>o depends on the boundary magnetic field h (at h = we have 4>o = 0). These boundary 
conditions preserve the topological U (1) symmetry giving conservation of soliton number, 
but they break the Z2 charge conjugation symmetry relating soliton and antisoliton (unless 
h = 0). 

Because we are mainly interested in checking the physical S matrices and paving 
the way towards a thermodynamic analysis, we shall discuss the lattice model and the 
underlying inverse scattering problem in the presence of boundaries as little as possible. 



We refer the reader to fll5| , |TT|| for details. We shall simply use the resulting bare Bethe 
ansatz equations without deriving them. 
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2. The bare Bethe ansatz equations 

The wave function of the inhomogeneous six- vertex model can be expressed in terms 
of a set of "roots" ctj, where j = 1 . . . n. They must be solutions of the set of equations 

N [f(aj + A, 7) + f(aj - A, 7 )] + 2f(a hl H) = 



where lj is an integer. The function / is defined as 

smh{ib — a)/2 



(3) 



f(a, b) = —i In 
and 



sinh(i6 + a)/2 



H = -fii'y, -iln(/i + cos7)). (4) 

7 

For h = 0, 7i? = 7r — 7. By construction of the Bethe-ansatz wave function, aj > 0. Even 
though there is a solution of (f|) with one vanishing root for any N and n, we emphasize 
that ctj = is not allowed because the wave function vanishes identically in this case. 

The solutions of these equations are quite intricate for arbitrary 7 [fL8t . For simplicity, 
we restrict to the case 7 = n/t where t is an integer, and consider both choices e = ±1. 
In the sine-Gordon model, these fall in the attractive and repulsive regimes respectively. 
We make the standard assumption that all the solutions of interest are collections of u k- 
strings" for k = 1, 2 . . A— 1 and antistrings a jn|. A /c-string is a group of aj in the pattern 

— in(k — 1), — irc(k — 3), ... , a' k ' +in(k — 1) where a" is real. The antistring has 
otj = + Z7r, where aW is real. We shall comment on the validity of this assumption 
(which is more crucial when one considers subleading effects) in the conclusion. 

The thermodynamic limit is obtained by sending A^ — > 00. In this case, we can define 
densities of the different kinds of solutions. The number of allowed solutions of (0) of type 
k in the interval (a, a + da) is (pfe(a) + p^(a))2Nda, where pk is the density of "filled" 
solutions (those which appear in the sum in the right-hand-side of (§) ) and p\ is the 
density of "holes" (unfilled solutions). The densities p a and p 1 ^ are defined likewise for the 
antistring. The "bare" Bethe ansatz equations follow from taking the derivative of (|3|). 
For 7 = n/t they can be written in the form: 

1 l ~ l 1 

2n(p k + pi) = - [o fc (a + A) + a k (a - A)] - (f) k ,t-i * p a + V" </>/tz * pi + 7^r T u k 

1 1=1 ^ (5) 

Mpa + Pa) = Mpt-l + Pt-l) + ^{U a ~ Ut-l) 



where * denotes convolution: 



/oo 
da'f(a - a')g(a') 
-oo 



These densities are originally defined for a > 0, but the equations allow us to define 
Pk(— ot) = Pk((x) in order to rewrite the integrals to go from — oo to oo. The kernels in 
these equations are defined most easily in terms of their Fourier transforms: 



oo 



I da infx /~ , , \ . cosh a; smh(t — k)x smhZx 

<f) kl (x) = / — e iatx/7T <f)ki(a) = S ab - 2 - 1 , (6) 

J_ oq 2tv sinh x sinh tx 

for k > I with <pik = tfiki, and 

sinh(t — k)x 



a-k = 



sinh tx 

sinh(t — H)x sinh kx sinh(t — 2k)x/2 



k sinh a; sinh sinhtx/2 ^ 

sinhifx sinh(t — 2)x/2 

u a = 2 — ^ — y- 1. 

sinhtx sinhtx/2 

The boundary manifests itself in the first term in Uk ; notice that even for h = it still 
modifies the equations. A few technicalities account for the other terms (these are relevant 
here because we are interested in subleading boundary effects). The second term in Uk 
arises from the fact that the sum in (^|) does not include the term m = j; the integration 
over densities includes such a contribution and so it must be subtracted off by hand. The 
third term in Uk arises because p and p h are defined for allowed solutions, while as already 
explained, a = is not allowed because it does not give a valid wavefunction. Since it is 
a valid solution of (|^) but is not included in the densities, we must subtract an explicit 
^fj5(a) (corresponding to 1/2N in Fourier space). One expects all other corrections to be 
suppressed by additional powers of 1/N. 

3. The physical Bethe ansatz equations 

We now derive the physical Bethe ansatz equations by extending the computation 



done in |T^] for periodic boundary conditions. The two cases e = ±1 are of course very 
different. 

Case i: e = 1 

In this case the ground state is filled with anti-strings so that pk = and p% = 0. We 
can use (0) to find the densities in the ground state, but we will not need these explicitly 



here. Excited states are given by including string solutions (p k 7^ 0) and/or holes in the 
antistring distribution (pjj 7^ 0). The physical equations are obtained by rewriting the 
bare equations ((5|) so that physical densities (i.e. p k and p 1 ^) appear on the right-hand side. 
This is done simply in Fourier space by solving for the Fourier transform p a in the last 
equation in (H) and substituting it into the others. The result is 

1 t_1 1 

27v(p k + p h k ) = - [A k {a + A) + A k (a - A)] + * p h a + $w * Pi + ^U k 

1=1 (8) 

where the kernels this time are 

~ . coshx coshft — 1 — k)x sinh/x , , . , 

$kl = hi - 2 \ .' k, I ^ t - 1; k > I 

cosh(r — l)x sinhx 

^ coshx sinh /cx ^ _^ ^ ^ 

cosh(t — l)x sinhx 1 

~ sinh(£ — 2)x 

$t-i,t-i ^ 



2 cosh(t — l)x sinhx 
where $^ = and. 

A k = cosh( f, : 1 ~ k)x k*t-i 

cosh(r — ljx 

A - 1 
11 2 cosh(t — l)x 

cosh(t — 1 — H)x sinh kx cosh x sinh kx /sinh(t — 2)x/2 

U k = 2 — I -(- 1 

cosh(t — l)x sinhx cosh(t — l)x sinhx \ sinhtx/2 

sinhft - 2/c)x/2 , , 

H H ^r 1 1 k^t-1 

sinh tx/2 

sinh(2t - 2 - ff)x sinh(t - 2)x/2 cosh fx sinh(t - 2)x 



(10) 



cosh(t — l)x sinhx cosh(t — l)x sinhx 2 cosh(t — l)x sinhx 
^ sinhifx sinh(t — 2)x/2 cosh f x sinh(t — 2)x 

a cosh(i — l)x sinhx cosh(t — l)x sinhx 2 cosh(t — l)x sinhx 

(We did not simplify some of the trigonometric sums to allow further identification of 
the various terms.) These are the physical Bethe ansatz equations, governing how the 
actual quasiparticle excitations interact with each other. Notice that the kernels are now 
symmetric, as opposed to ([5]). Each density corresponds to a quasiparticle; it is easy to 
read off the bulk S matrix elements. In the next section, we will discuss how to do this 
and how to find the boundary S matrix. 
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Case %%: e = — 1 

In this case the ground state is made of real solutions. The physical densities are 
therefore p\, Pk and p a . Eliminating p\ from the right hand side of the bare equations 
gives 

2ir( Pl + p\) = hs(a - A) + s(a + A)] + s * of _1) * p\ 



ajfc_i * Pk + a\_ 2 *p a + —u 



27T(p, + p%) = 4_ x * + ^ j_! * - 4-l't-l * Pa + 2N U k-l 

1=2 



MPa + Pa) = MPt-l + Pt-l) + ^ 



where the quantities ^ and u^_}' H ^ are the same as in the bare equation 

®"dZD whh t,k,H there replaced by £ — l,k — 1, H — 1 respectively, and 

sinhft — H)x coshtx/2 sinhft — 2)x/2 sinhft — 2)x 

u = ; : T71 7\ H i T\ H 



coshx sinh(t — l)x cosh x sinh(t — l)x 2 coshx sinh(t — l)x 

1 



2 cosh x 

It is more difficult to read off the quasiparticle spectrum and the bulk and boundary S 
matrices here because the bulk scattering is non-diagonal and most densities correspond 
in fact to pseudoparticles. We will discuss this in the next section. 

4. The S Matrix 

In this section we derive the physical Bethe ansatz equations from the bulk and bound- 
ary S matrices by quantizing the momenta of a set of relativistic quasiparticles on a line. 
We identify these with the physical equations derived from the lattice model in the cases 
e = ±1. 

In a relativistic quantum field theory the energy and momentum of a particle of mass 
m can be parametrized in terms of the rapidity 9: (E, P) = (m cosh 9, msinh#). Lorentz 
invariance requires that S matrix elements depend only on rapidity differences. 

First let us discuss the case when the S matrix is diagonal. When there are p different 
species of particle, we have the two-particle S-matrix elements Sb c {9\ — 6> 2 ), where b and 
c run from 1 to p. These S'-matrix elements give the phase shift in the wavefunction 
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when two particles are exchanged. We also have the boundary S matrix elements R b {9) 1 
which gives the phase shift when a particle of species b bounces off a wall and changes its 
rapidity from 9 to —9. (We do not consider the case where a particle changes species when 
bouncing off the boundary.) We have a gas of M particles on a line of length L, with the 
ith particle of species 6j. For the fixed boundary conditions we consider, we demand that 
the wavefunction vanish at both ends of the line. This requires making a stationary wave 
of states with opposite momenta, with the momenta subject to the constraints 

M 

e *Lm 6i sinhfl, "Q S bibj (9 i -9 j )R bi (9 i ) 

\\ Sbibj (~0i - 0j)R bi (-9i) 



—iLnrifr . sinh 6f 

— C 1 



for all i. When the scattering is trivial, this reduces to the familiar P = mr/L, where n is 
a positive integer. Unitarity requires that S~ 1 (9) = S(—9) and R~ 1 (9) = R(—9), so this 
can be rewritten as 



^2im b . sinh 9iL 



J] S hbj (9 l -9 J )S bibj (9 l + 9 J )Rl(9 l ) = l. (13) 

Alternatively we can deduce these equations from systematic application of the Zamolod- 
chikov-Faddeev algebra together with its boundary counterpart. 

As in the lattice model, we define densities of solutions of (|IB|) so that (p b (9) + 
p b l (9))2Ld9 gives the number of allowed rapidities between 9 and 9 + d9 for species b, 
and p b gives the density of filled states. Taking the logarithm of (|T3| ) gives an equation of 
the form (|||); taking the derivative of this gives an equation for the densities. We again 
define the densities at negative rapidity by p b (—9) = p b {9). We have one equation for 
every type of particle: 

27r(p b + p^) = m b cosh 9 + ^ Pbc * Pc + 7Ty®b, (14) 



where 



2L 

c=l 



<Pbc(0) = -i-^g lnS bc (9) 

e 6 (0) = -i^-2\nR b (9) + i^-\nS bb (29) - 2ivd(9) 
(19 d9 



(15) 



The additional terms in 0;, result from the same considerations as in the lattice model: 
the facts that the product in ( |13D has no term i = j and that the solution with 9 = 
should not be included. As before, we expect additional corrections to be suppressed by 
powers of 1/L. 

The equations ([14]) are the physical Bethe ansatz equations derived directly from 
the S matrix description. They enable us to compare a diagonal S matrix to lattice 
results. We do this in case i below. If the scattering is not diagonal, the analysis is more 
complicated: one must "diagonalize" the multiparticle state, including the boundary S 
matrices. Basically, one must find the eigenvalues A(9{\{9j}) of the object 

T(9 l \{9 J }) bc R c (9)T- 1 (-9 l \{9 J }) cb R^ 1 (-9 l ) (16) 

where T(9i\{9j})b c is the matrix which describes scattering a particle of species b and 
rapidity 9i through all the other particles and ending up with that particle being of species 



c. The generalization of (|13|) then is 



e 



limy, . sinh 9i L 



m\{0j}) = i. (17) 



To diagonalize ([LB)), one uses the Bethe ansatz and the analysis proceeds as before, but the 



calculation is a little delicate. It is however simplified by using results of Hlq] ; we discuss 
this in case ii below. 

Case i: f3 2 = 8n/t 

As discussed in the introduction, the scaling limit of the inhomogeneous XXZ model 
with e = 1 is the sine-Gordon model at /3 2 = 8n/t. The particles and their bulk S 
matrix in this case are well known ]l!| . There are t — 2 breather states with mass m a = 
2M sinh 2 {t-i) ' a l° n g with a soliton and an antisoliton of mass M. All of these states scatter 
diagonally with each other. The boundary S matrices for the soliton and antisoliton are 
conjectured in 0], and those for the breathers in 0. 

Using these S matrix elements, we can now compare the physical Bethe ansatz equa- 
tions (0) with the lattice ones (§) . First we have to extract the scaling limit of the lattice 
model. This is easily done by taking the limit A — ► oo at finite bare rapidity a. The source 
term in equations (H) reproduces the expected spectrum with 

M OC l e -*A/2(t-l) 

A 



and the physical rapidity 
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-a, 



2(t-l) 

provided we identify the breathers with the string solutions k = 1 . . . t — 2, the soliton with 
the string solution t — 1, and the antisoliton with the antistring a. 

It is already known JT2J that the kernels $jfc in (|9|) are the same as the p^c that result 
from the known Sb c scattering matrix. It is straightforward but rather tedious to check 



that the kernels Uk in ( pTOl) too are the same as ©t, resulting from ([15]) and the Sb c (9) and 
Rb(9) of ||T9|,3,|6| , completing the boundary identification. Let us discuss the case of soliton 
and antisoliton only for simplicity, denoting them + and — respectively. Using the integral 
representation for the log of a gamma function, the result of jl can be written 



-i\nR±{6) 



dx . 2 

— sin — (t — l)xu 

x 7T 

sinh(t - 1 ± 2£ > /n)x sinh |(t - l)xsinh(| - l)x 
2cosh(i — l)xsinha; sinhx/2 sinh2(t — l)x 



where £ depends on 



io- Comparing this with (|T0| ) allows us to identify 

To identify the remaining terms we recall that 

sinh(t 



-ilnS++(0) 



dx sinhft — 2)x 2 , 

K —- — sin -it - l)xO. 

x 2 smhx cosh(t — l)x 7r 



The identity 



sinh3(t - l)x/2sinh(t - 2)x/2 sinh(t - 2)x/2coshte/2 



sinhx/2 sinh2(t — l)x 

sinh(t - 2)x/2 



2 cosh(i — l)x sinh x 
sinh(i — 2)x 



(18) 



+ 



4 sinhx/2 cosh(t — l)x/2 4 sinh a; cosh(t — l)x 

completes the desired identification. 

This confirms the £ matrices conjectured in at these values of (3. We emphasize 
that it is absolutely crucial to have included the correction terms to Uk in (^) and in 
( |l~5D to get this result. 

Case ii: ft 2 = 8n(t - l)/t 

The continuum limit of the lattice model at 7 = Ti/t is the sine-Gordon model at 
(3 2 = 8n(t — l)/t. The only states in the spectrum are the soliton and antisoliton, and they 
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scatter with a nondiagonal S matrix [19]. With our fixed boundary conditions, soliton 
number is conserved, so the scattering off the wall is diagonal but only if h = (H = t — 1) 
is it equal for soliton and antisoliton ||. 

To derive the Bethe ansatz equations for this nondiagonal S matrix, we need the eigen- 
values of the matrix ([16]) . We now discuss some technical aspects of this computation; the 
uninterested reader can jump directly to the next paragraph. The equivalent problem for 
periodic boundary conditions, where one diagonalizes T^b, is well understood. In this case, 
a simplifying feature is that S(0) is simply the permutation operator. Hence, instead of Tbb 
one can consider the matrix describing the scattering of one particle through all the others 
including itself. The unwanted contribution of a particle scattering with itself is trivial 
anyway and does not have to be discarded. The advantage is that the new diagonalization 
problem fits perfectly in the inverse scattering framework and is easy to solve. In the 
presence of a boundary the situation is different. If one tries the same approach there are 
now two unwanted contributions, since the particle at rapidity 9 can "scatter with itself" 

at rapidity ±9. The process with S(9 9) = S{29) leads to non-trivial terms. When 

the scattering was diagonal this was easy to discard by hand, but the situation is more 
complicated here. On the other hand, the inverse scattering formalism in the presence of 
boundaries developed e.g. in ]15|,|17| * s treating a slightly different problem: R b ~ 1 (—9) is 
replaced there with R^ 1 {—9 — in). Fortunately these two complications essentially cancel 
each other; changing the rapidity — 9 to —9 — in acts as a projector which cancels out 
the unwanted processes resulting from the self-scattering term. The final expression for 
the eigenvalue A(9i\{9j}) turns out to be the same as that of |y|, up to some missing 
prefactors that produce a phase, as required. The complete proof is rather tedious and 
will be provided upon request. The interested reader can also check the final result by 
direct calculation in the two-particle case. 

The result is that the auxiliary problem determining the eigenvalue A is similar to the 
bare problem, up to a redefinition of some parameters. Let us define u{9), which is the 
density of particles (soliton or antisoliton). In this case we have 



-/•ln/?i(0)= / — 

x 



sinh(l ± 2£(t - 1) /tt)x sinh \x sinh(t - 2)x/2 
2 coshx sinh(i — l)x sinh(t — l)x/2 sinh2x 



2 

sin — x9 

7T 



together with 



dx sinh(t-2)x . 2 

-% mb ++ (9) = — / : — — sm 

x 2 coshx smh(t — lja; tv 
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The first part of the quantization equation (the derivative of the log of (|13|)) reads 

i d 

2r:(a + a h ) = mcosh#- — — In A, (19) 

2L do 

with 

" 5ll ln A = lnS++(9) ) * " - 2 of" 1 ' * « + J& * a. + ^, 

7 b=l 



and 



17 = hi !?+(#) + ^ l n> S ++ (2^) - 2n5(8). 

do do 



For the other densities aj, and (b = 1 . . . t— 2), we have the bare equations @ with three 
changes. First, t is replaced by t — 1 because there are t — 1 solutions to the new Bethe 
equations. There is also a new source term 1 * a because the "vertices" in the problem 



of diagonalizing (|16|) correspond to the real particles, and have fluctuating numbers and 
rapidities. These two changes are identical to the case with periodic boundary conditions. 



In addition, we have a boundary term which from JT^] is the same as the one in the bare 
equations with H replaced by H — 1. This shift is of course related to the shift of t. For 
instance, if in the original bare equations h = so that H = t — 1, then in the auxiliary 
problem we again expect soliton-antisoliton symmetry yielding H = (t — 1) — 1. 
We can thus identify ( |19D and the equations for Ob with (|IT]) provided we set 

7T / H 



2 V t-1 

by using again the identity (|T8| ) with x replaced by (t — l)x and t — 1 by l/(t — 1). This 
confirms results of 0|| at values of (3 where the scattering is non-diagonal. 

5. Conclusion 

We have verified in both attractive and repulsive regimes that the boundary S matrices 
conjectured in [Q,|6|] are in fact correct. We think that the main interest of our computation 
lies in the delicate analysis of all the boundary terms, which are crucial to make the 
computation consistent. In particular, observe from the example of a = that at order 
l/N, the modification of even a single root in the Bethe equations affects the results. This 
puts the string hypothesis to a much stronger test than the bulk computations done so far. 
If we believe that the picture of 0] is consistent and that the inhomogeneous XXZ model 
provides also a regularization of the sine-Gordon model in the presence of boundaries, we 
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are forced to conclude that strings are the only solutions of the Bethe equations that 
contribute to the scaling limit, i.e. to the low-energy excitations. This does not contradict 



the known examples of non-string solutions, which either exist in finite numbers [|10[|20|1 
(so we show that the finite number is zero in the scaling limit) or are at very high energy 

The next step is to study the system at finite temperature and to derive the thermo- 
dynamic Bethe ansatz equations for the free energy. This should yield the "(/-function" , 
which measures the number of degrees of freedom on the boundary ||22|| . At a fixed point, 
this can be calculated from conformal field theory p3|,p2|, and was calculated for the free 



boson (the critical limit of the sine-Gordon model) in [24|]. It is simple to get the bulk 
part of the free energy from the physical Bethe ansatz equations, but additional subtleties 
appear because the g-function is a subleading term. For example, the usual thermody- 
namic Bethe ansatz free energy is the log of the partition function only at leading order 
in L; there are system-dependent corrections at lower order. However, when one looks at 
massless S matrices to describe a flow between ultraviolet and infrared fixed points, the 
ratio guv / 9m should be independent of these corrections. It is indeed simple to calculate 
this ratio from the S matrix, for example, in the Kondo problem 0,^5] or in the flows 
between the minimal models. We hope to discuss these more physical aspects in the near 
future. 



Acknowledgments: This work was supported by the Packard Foundation, the 
National Young Investigator program (NSF-PHY-9357207) and the DOE (DE-FG03- 
84ER40168). We thank A. Delfino and G. Mussardo for very useful discussions on boundary 
problems. 



12 



References 



[1] For a review, see I. Affleck, "Conformal Field Theory Approach to Quantum Impurity 

Problems", UBCTP-93-25, |cond-mat/9311054 . 
[2] N. Andrei, K. Furuya, and J. Lowenstein, Rev. Mod. Phys. 55 (1983) 331; 

A.M. Tsvelick and PB. Wiegmann, Adv. Phys. 32 (1983) 453. 
[3] I. Cherednik, Theor. Math. Phys. 61 (1984) 977. 

[4] S. Ghoshal and A.B. Zamolodchikov, "Boundary State and Boundary S Matrix in 
Two-Dimensional Integrable Field Theory", RU-93-20, |hep-th/ 9306002. 



[5] A. Fring and R. Koberle, "Factorized Scattering in the Presence of Reflecting Bound- 
aries", USP-IFQSC/TH/93-06, |hep-th/9304141 . 



[6] S. Ghoshal, "Bound State Boundary S Matrix of the Sine-Gordon Model", RU-93-51, 
[hep-th/9310188 . 



[7] A. Fring and R. Koberle, "Affine Toda Field Theory in the Presence of Reflecting 

Boundaries", USP-IFQSC/TH/93-12, |hep-th/9309142| . 
[8] R. Sasaki, "Reflection Bootstrap Equations for Toda Field Theory" , [hep-th/9311027" . 
[9] S. Ghoshal, "Boundary S Matrix of the 0(n) Symmetric Nonlinear Sigma Model" 
RU-94-02, frep-th/94010Qg . 
[10] C. Destri and J.H. Lowenstein, Nucl. Phys. B205 (1982) 369. 
[11] G.I. Dzapardize, A. A. Nersesyan and P.B. Wiegmann, Phys. Scr. 27 (1983) 5. 
[12] A.N. Kirillov and N. Reshetikhin, J. Phys. A20 (1987) 1565, 1587. 
[13] C. Destri and H. de Vega, J. Phys. A22 (1989) 1329. 

[14] N. Reshetikhin and H. Saleur, "Lattice Regularization of Massive and Massless Field 
Theories", USC-93-020, |hep-th/9309i35 . 



[15] E.K. Sklyanin, J. Phys. A21 (1988) 2375. 

[16] F. Alcaraz, M. Barber, M. Batchelor, R. Baxter and G. Quispel, J. Phys. A20 (1987) 
6397. 

[17] L. Mezincescu and R.I. Nepomechie, Int. J. Mod. Phys. A6 (1991) 5231; Int. J. Mod. 

Phys. A7 (1992) 565; H.J. de Vega and A. Gonzalez-Ruiz, preprint LPTHE-93/38. 
[18] M. Takahashi and M. Suzuki, Prog. Th. Phys. 48 (1972) 2187. 
[19] A.B. Zamolodchikov and ALB. Zamolodchikov, Ann. Phys. 120 (1980) 253. 
[20] F. Woynarovich, J. Phys. A15 (1982) 2985; 

O. Babelon, H.J. de Vega and CM. Viallet, Nucl. Phys. B220 (1983) 13. 
[21] F. Essler, V. Korepin and K. Schoutens, J. Phys A25 (1992) 4115. 
[22] I. Affleck and A. Ludwig, Phys. Rev. Lett. 67 (1991) 161. 
[23] J. Cardy, Nucl. Phys. B324 (1989) 581. 
[24] S. Eggert and I. Affleck, Phys. Rev. B46 (1992) 10866. 
[25] P. Fendley, Phys. Rev. Lett. 71 (1993) 2485, |cond-mat/9304031 . 



13 



